MINIMAL CUBATURE RULES ON AN UNBOUNDED DOMAIN 
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Abstract. A family of minimal cubature rules is established on an unbounded 
domain, which is the first such family known on unbounded domains. The 
nodes of such cubature rules are common zeros of certain orthogonal polyno- 
mials on the unbounded domain, which are also constructed. 



1. Introduction 

In two or more variables, few families of minimal cubature rules are known in 
the literature, none on unbounded domains. The purpose of this note is to record 
a family of minimal cubature rules on an unbounded domain. 

The precision of a cubature rule is usually measured by the degrees of polynomials 
that can be evaluated exactly. For a nonnegative integer m, we denote be n 2 ^ the 
space of polynomials of degree at most m. Let fi be a domain in R 2 and let W 
be a non- negative weight function on f2. A cubature rule of precision 2n — 1 with 
respect to W is a finite sum that satisfies 

N 

(1.1) / f(x,y)W(x,y)dxdy = y2\ k f(x k> y k ), V/en^, 

Jn fc =i 

and there exists at least one function /* in n|„ such that the equation (1.1 1 docs 
not hold. A cubature rule in the form of (JTTTJ) is called minimal, if its number of 
nodes is the smallest among all cubature rules of the same precision for the same 
integral. 



It is well known that the number of nodes, N, of (1.1) satisfies (cf. [HE]), 

(1-2) iv>dimnL^^±ll. 

A cubature rule of degree In — 1 with N = dimR 2 ^ is called Gaussian. In contrast 
to the Gaussian quadrature of one variable, Gaussian cubature rules rarely exists 
and there are two family of examples known [31 [7], both on bounded domains. It 
is known that they do not exist if W is centrally symmetric, which means that 
W(x) = W(— x) and —x G f2 whenever x G f2. In fact, in the centrally symmetric 
case, the number of nodes of ( |1.1[ ) satisfies a better lower bound 0], 

(1.3) 7V>dimII 2 ■ n -< n+l ) 
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A cubature rule whose number of nodes attains a known lower bound is nec- 
essarily minimal. It turns out that a family of weight functions W Q p±x defined 

by 

(1.4) W a . JSt±i (x,y)^\x + y\ 2a + 1 \x~y\^+\l-x 2 ) ± Hl~y 2 ) ± K a,[3>-l, 

on [— 1, l] 2 admits minimal cubature rules of degree An — 1 on the square [— 1, 1] , 
which was established in [5] for a = (3 = (see also [TJ [TD] ) and in [TT] for (a, (3) ^ 
(0, 0). There are not many other cases for which minimal cubature rules are known 
to exist for all n and none in the literature that are known for unbounded domains. 

To establish the minimal cubature rules for W a> p ±12, the starting point in [11] 
is the Gaussian cubature rules in [7] and it amounts to a series of changing of 
variables from the product Jacobi weight function on the square [— 1, l] 2 to the 
weight function W a ,p,±i2- The procedure works for general product weight function 
on the square. Moreover, as we shall shown in this note, that the procedure also 
works for an unbounded domain, which leads to our main results in this note. The 
minimal cubature rules are known to be closely related to orthogonal polynomials, 
as their nodes are necessarily zeros of certain orthogonal polynomials. We will 
discuss this connection and construct an explicit orthogonal basis on our unbounded 
domain in Section 3. 

2. Gaussian Cubature rules on an unbounded domain 

Let W be a nonncgative weight function on a domain O C R 2 . A polynomial 
P G IT^ is an orthogonal polynomial of degree n with respect to W if 

/ Pk,n( x :y)<l( x i y)W(x,y)dxdy = 0, Vg G 
Jn 

Let V 2 be the space of orthogonal polynomials of degree exactly n. Then dim V 2 = 
n + 1. The Gaussian cubature rules can be characterized in terms of the common 
zeros of elements in V 2 ([5H5]). 

Theorem 2.1. Let {Pk, n '■ < k < n} be a basis ofV„. Then a Gaussian cubature 
rule of degree 2n — 1 for the integral against W exists if and only if its nodes are 
the common zeros of Pk, n > < k < n. 

We now describe a family of Gaussian cubature rules on an unbounded do- 
main. Let w(x) be a nonnegative weight function defined on the unbounded domain 
[1, 00) and let c w denote its normalization constant defined by c w J ± w{x)dx = 1. 
Let p n (w;x) be the orthogonal polynomial of degree n with respect to w and let 
X2,m ■ ■ ■ , x n _ n be the zeros of p n (w] x). It is well known that Xk, n are real and 
distinct points in [l,oo). The Gaussian quadrature rule for the integral against w 
is given by 

/oo n 
f(x)w(x)dx = ^ Afc,n/(£fe,n)) / G n 2n -l, 
fe=l 

where n 2 „_i denotes the space of polynomials of degree 2n — 1 in one variable and 
the weights \k,n are known to be all positive. 

A typical example of w is the shifted Laguerre weight function 

w a (x) := (x-l) a e~ x+1 , a>-l, 
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for which the orthogonal polynomial p n (w a ; x) is the Laguerre polynomial L"(x— 1) 
with argument x — 1. 

The unbounded domain on which our Gaussian cubature rules live is given by 

(2.2) n := {(u,v) : < u- 1 < v < u 2 /A}, 

which is bounded by a parabola and a line for u > 1 and v > 2, which is the shaded 
area depicted in Figure 1. 




Figure 1. Domain £1 



The function w(x)w(y) is evidently a symmetric function in x and y. For 7 > — 1, 
we define the weight 

(2.3) W y (u,v) := w(x)w(y)\u 2 -4«| 7 , 7 > -1, 
where the variables (x,y) and (u, v) are related by 

(2.4) u = x + y, v = xy. 

The function w(x)w(y) is obviously symmetric in x and y, so that it can be written 
as a function of (it, w) for (x, y) in the domain 

A := {(x,y) : 1 < x < y < 00}. 

The function W~(u,v) is the image of w{x)w{y)\x — y| 27+1 under the changing of 
variables u = x + y and v = xy, which has a Jacobian |x — y\ and |x — y| = \J v? — 4v. 
When w = w a is the shifted Laguerre weight, we denote W 1 by W an , which is given 
explicitly by 

W an (x,y) := (v-u + l) a e- u - 2 (u 2 - 4w) 7 . 
The changing of variables (|2.4|) immediately leads to the relation 



(2.5) / f(u,v)W 1 (u,v)dudv= / f(x + y, xy)w(x)w(y)\x — y\ 2l+1 dxdy 



A 



[1,00) 



f(x + y, xy)w(x)w(y)\x - y\ 2l+1 dxdy, 
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where the second equation follows from symmetry. Now, recall that x k , n denote 
the zeros of p n (w;x). We define 

^k,j — X k _ n ~l~ Xj : ji 7 ^k,j — X k ^ n Xj^ nj — J — k — ^l- 

Theorem 2.2. For W_i on Q, the Gaussian cubature rule of degree 2n — 1 is 
(2.6) c 2 w I f(u,v)W_i('u,v)dudv = ^2^2 X k Xjf(u jtk ,v jik ), f e n| n _ 1; 

•' n k=l j=l 

where ^ means that the term for j = k is divided by 2. For Wi on Q, the Gaussian 
cubature rule of degree In — 3 is 

„ n fc— 1 

( 2 -7) c 2 w f(u,v)Wi(u,v)dudv = ^2^2x jtk f(uj,k,v j>k ), /en^ 3 , 

Jn ' k=2j=l 

where X jjk = XjX k (xj in - x k . n ) 2 . 

The proof follows almost verbatim from the proof of Theorem 3.1 in [11] , In fact, 



by (2.5), the cubature rule (2.6) is equivalent to the cubature rule for w{x)w{y)dxdy 
on [l,oo) 2 for polynomials f[x + y,xy) with / £ Il| ra _ 1 , which are symmetric 
polynomials in x and y. By the Sobolev's theorem on invariant cubature rules [5], 
this cubature rule is equivalent to the product Gaussian cubature rules for w(x)w(y) 



on [1, oo) 2 , which is the product of (2.1 ). Thus, the cubature rule (2.6) follows from 



the product of Gaussian cubature rules under (2.4). 

The existence of these cubature rules also follow from counting common zeros of 
the orthogonal polynomials with respect to W±i. Indeed, a mutually orthogonal 
basis with respect to W_i on Q is given by 

(2.8) P k \ n k) (u,v)= Pn {x)p k (y)+p n (y)p k (x), < k < n, 
and a mutually orthogonal basis with respect to Wi is given by 

(2.9) p (h) {UiV) = Pn + i{x)Pk(y)- Pn+ MPk{x) ^ < fc < n; 

x — y 



both families are defined under the mapping (2.4). This was established in [2] for 



the domain [—1,1] , but the proof can be easily extended to our Q,. It is easy to 
see that the elements of {(%,&, Vj,k) '■ < j < k < n} are common zeros of P k n 2 , 
< k < n, and the cardinality of this set is dimll 2 ^, which implies, by Theorem 
2.1, that the Gaussian cubature rules for W_i exists. The proof for Wi works 

2 2 

similarly. 

3. Minimal cubature rules on an unbounded domain 
We are looking for cubature rules of degree 2n — 1 that satisfy the lower bound 



(1.3), which are necessarily minimal cubature rules. Such cubature rules are char- 



acterized by common zeros of a subspace of V 2 ( [1] ) • 

Theorem 3.1. A cubature rule whose number of nodes attains the lower bound 
(1.3) exists if, and only if its nodes are common zeros of L^^J + 1 orthogonal 



polynomials of degree n. 
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Let w be the weight function defined on [1, oo) and let W 7 be the corresponding 
weight function in (2.3) defined on £1 in (2.2). We define a family of new weight 
functions by 

(3.1) W 7 (x,y) ~W 7 (2xy,x 2 +y 2 -l)\x 2 -y 2 \, (x,y)eG, 

where the domain G is centrally symmetric and defined by 

G := (-oo,-l] 2 U [l,oo) 2 . 



Thus, the weight function W 7 is centrally symmetric on G 

o 

A = {(x, y) : 1 < x < y < oo} 



That W 7 is well-defined on G is established in the next lemma. Recall that 



Lemma 3.2. The mapping (x, y) i— > (2xy, x 2 + y 2 — 1) is a bijection from A onto 
n. Furthermore, 

(3.2) / f(u,v)W 7 (u,v)dudv = / f{2xy,x 2 + y 2 -l)W 1 (x,y)dxdy. 
Jn jg 

Proof. Recall that if (x, y) € A, then (x+y, xy) G and the mapping is one-to-one. 
For (x, y) G [1, oo) 2 , let us write x = cosh# and y = cosh<^>, < 6, <p < n. Then it 
is easy to verify that 

(3.3) 2xy = cosh(6»-0)+cosh(6' + 0), x 2 + y 2 - 1 = cosh(6< - cf>) cosh(6l + 0), 

from which it follows readily that (2xy, x 2 + y 2 — 1) E ft whenever (x, y) G A. The 
Jacobian of the change of variables u — 2xy and v — x 2 + y 2 — 1 is A\x 2 — y 2 \, so 
that the mapping is a bijection. Since dudv — A\x 2 — y 2 \dxdy and the area of G is 
four times of A, the formula (3.3) follows from the change of variables, the integral 
(2.5) and the fact that f(2xy, x 2 + y 2 — 1) is central symmetric on G. □ 

In the case of W 7 = W a ,j, we denote the weight function W 7 by W a , 7 , which is 
given explicitly by 

(3.4) W a , 7 (x, y) := Ai\x - y\ 2a (x 2 - l)^(y 2 - ir\x 2 - y 2 ^ 2 ^ 2 . 

For the weight function W 7 on the unbounded domain G, minimal cubature rules 
of degree An — 1 exist, as shown in the next theorem. To state the theorem, we 
need a notation. For n G No, let Xk, n be the zeros of the orthogonal polynomial 
p n (w;x), which are in the support set [l,oo) of the weight function w. We define 
Ok.n by 

Xk.n = cosh9k, n , 1 < k < n. 
Since Xk n > 1, it is evident that Ok n are well defined. We then define 



(3.5) Sj.k ■= cosh 3 '" 2 k,n and tj t k :— cosh 



2 "},* • 2 

Theorem 3.3. For W_ i on G, we Ziawe t/ie minimal cubature rule of degree 4n— 1 
witfi dimll2 n _ 1 + n nodes, 

(3-6) c 2 ,, / f(x,y)W_i(x,y)dxdy =^5ZX1 ^^k t n[f(sj,k,tj,k)+f(tj,k,Sj,k) 

Jg k=l j=l 



+ f(Sj,k,-tj,k) + f(~tj,k, Sj,k)} ■ 
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ForWi on G, we have the minimal cubature rule of degree An— 3 with dimll 2 ,,^ 
nodes, 

(3.7) 



n k— 1 



/ y)Wi (x, 2/)dxefy =- 2J 2J A ^ if( s i,k,tj,k) + f(tj,k,Sj,k) 
Jg fe=2 i=i 

+ f(Sj,k,-tj,k) + f{-tj,k, Sj,k)] , 

where \ jt k = Aj,„Afe j „(cosli^, Ti - cosh6> fcj „) 2 . 

Proof. We prove only the case of W_ i , the case of Wi is similar. Our starting 
point is the Gaussian cubature rule in (2.6), which gives, by (3.2 1, 

k 



c w / f( 2x yi x2 +y 2 - 1 )^ } -^(x,y)dxdy = ^2^2 X j,n X k,nf(Uj,k,Vj,k) 
J G j — i — i 



By (3.3) or by direct verification 
cosh0j,„ + cosh^fe,„ : 

cosh 9j jTl cosh 9k, n = cosh 
which implies that 



k=l j = l 



2 cosh 'i."- 9 *-" cosh 9 ^ 6 



cosh 



2 flj,„+e fc , r . 



- 1 



3, n 



Xk,n 



and t>j ^ — Xj,n x k,n — s j.fe + t\k 



Consequently, the above cubature rule can be written as 



f(2xy, x 2 + y 2 - i (x, y)dxdy = J2Y, X k X jf( 2s iMi,k, Sj, k + tj, k - 1) 

' k=i i=i 

for all / € Il2„_ 1 . For / G Tl 2 n -i, tne polynomial f(2xyx 2 + y 2 — 1) is of degree 
4n — 1. Since the polynomials f(2xy,x 2 + y 2 — 1) are symmetric polynomials and 
all symmetric polynomials in Il| n _j can be written in this way, we have established 
(3.6) for symmetric polynomials. By the Sobolev's theorem on invariant cubature 
rules, this establishes (3.6) for all polynomials in H 2 n _ 1 . □ 



The number of nodes of the cubature rule in (3.6) is precisely 

2n 



AT = dim 



= dimll 



2n-l 



which is the lower bounded of (1.3) with n replaced by 2n. Thus, (3.6) attains the 
lower bound (1.3). Similarly, (3.7) attains the lower bound (1.3) with n replaced 
by 2n - 1. 

It turns out that a basis of orthogonal polynomials with respect to W 7 can be 
given explicitly. We need to define three more weight functions associated with W 7 . 

W\ hl) (u,v) := (l-u + v)(l + u + v)W~,(u,v), 

(3.8) W^'°\u,v) := (l-u + v)W^(u,v), (w,u)gO, 

W\ ' 1] {u,v) := (l + u + v)W 7 (u,v). 

Under the change of variables u = x + y and v = xy, 1 — u + v = (x — l)(y — 1) 
and 1 + u + v = (x + l)(y + 1). The three weight functions in (3.8) are evidently 



of the same type as W 1 . We denote by {Pj^ 



< k < n} an orthonormal 
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basis of V n (Wj) under 



P 



(7), 1,0, p(7),0,l 



'W-, 



For < fc < n, we further denote by P, 



(7), 1,1 



fe , ri- 



fe, n 



fc,n 



the orthonormal polynomials of degree n with respect to (/, g) 



li- 



fer w = wy , w 7 i,u; , w. 



(1,0) w (0,l) 



7 



respectively. 



Theorem 3.4. For n = 0, 1, . 



a mutually orthogonal basis o/V2„(H ; 7 ) is given 



,(7) 



(3.9) 



0< fc < 



and a mutually orthogonal basis of V2 n +i(VV 7 ) is given by 



1(7) 



(3.10) 



(7), 0,1 



(2xy,x 2 + y 2 -l), 



S„+i(z, 2/) := 0* " 2/)Pg' 1 '°(2a ; y, z 2 + y 2 - 1) 



< k < n- 1, 

< fc < n, 
< fc < n. 
(1.1)/ 



Under the mapping (u, v) H> (2xy, x + y — 1), it is easy to see that W-y ' (it, i>) 
becomes (x 2 — i/ 2 ) 2 W 7 (:k, y), Wj'°\u, v) becomes ( x— y ) 2 W-y(x, y), and W^ 0,1 \u, v) 
becomes [x + y) 2 W 1 (x, y). Hence, using Lemma 
lit 

similar to the proof of Theorem 3.4 in 



3.2 



the proof can be deduced 

from the orthogonality of Pk, n^' 1 ^ and symmetry of the integrals against w\ l, ^\ 



Combining with (2.8) and (2.9), we can express orthogonal polynomials iQ 
in terms of orthogonal polynomials in one variables. For example, in the case of 
Wa,7 in (3.4), we can express even degree orthogonal polynomials in terms of the 
Laguerre polynomials. 

Proposition 3.5. Let a > —1. A mutually orthogonal basis of V2n(W Q _ 1 ) is 
given by, for < fc < n and < k < n — 1, respectively, 

iQjQ^ ) (cosh0+ l,cosh0 + 1) = l4 Q) (cosh((9 - </>))4 Q) (cosh(6> + 0)) 

+ 4 Q) (cosh(0 - 0))4 a )(cosh(0 + 0)), 

2Q [ kfJ\ C0 ^ e + l,coBh^+ 1) = (x 2 - y 2 ) [4 a + 1) (cosh(^ - (j>))L [ * +1 \cosh{6 + 0)) 

+4 Q+1) (cosh(# - 0))4 t r f 1 1) (cosh(e + <j))) . 

In these formulas we used x = cosh# — 1 and y = cosh 4> — 1, so that L"(x — 
1) = i"(cosh#) and L"(y — 1) = L"(cosh</>) and we can then use (3.3). Note, 
however, we cannot express the odd degree ones in terms of Laguerre polynomials. 
In fact, for 2Q/c,2n+i, we need orthogonal polynomials of one variable with respect 
to w(x) = (x + i)(x — l) Q e~' :c_1 \ which is not a shift of the Laguerre weight 
function. 

Corollary 3.6. For the weight function W_i, the nodes of the minimal cubature 



formulas (3.6) are common zeros of orthogonal polynomials {iQ k 2 2 n '■ < k < n} 



An analogue result can be stated for the minimal cubature formulas (3.7). 
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